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Stellar collapse and the subsequent development of a core-collapse supernova explosion emit bursts 
of gravitational waves (GWs) that might be detected by the advanced generation of laser interfer- 
ometer gravitational-wave observatories such as Advanced LIGO, Advanced Virgo, and LCGT. 
GW bursts from core-collapse supernovae encode information on the intricate multi-dimensional 
dynamics at work at the core of a dying massive star and may provide direct evidence for the yet 
uncertain mechanism driving supernovae in massive stars. Recent multi-dimensional simulations 
of core-collapse supernovae exploding via the neutrino, magnetorotational, and acoustic explosion 
mechanisms have predicted GW signals which have distinct structure in both the time and frequency 
domains. Motivated by this, we describe a promising method for determining the most likely explo- 
sion mechanism underlying a hypothetical GW signal, based on Principal Component Analysis and 
Bayesian model selection. Using simulated Advanced LIGO noise and assuming a single detector 
and linear waveform polarization for simplicity, we demonstrate that our method can distinguish 
magnetorotational explosions throughout the Milky Way {D < lOkpc) and explosions driven by 
the neutrino and acoustic mechanisms to D < 2kpc. Furthermore, we show that we can differenti- 
ate between models for rotating accretion-induced collapse of massive white dwarfs and models of 
rotating iron core collapse with high reliability out to several kpc. 

PACS numbers: 04.30.Tv, 04.80.Nn, 05.45. Tp, 97.60. Bw 



I. INTRODUCTION 

Almost eighty years after the proposal by Baade & 
Zwicky that (core-collapse) supernovae represent the 
transition of an ordinary massive star into a neutron star 
[T], we still lack a complete understanding of this phe- 
nomenon. In particular, we do not know with certitude 
how the supernova mechanism operates and converts the 
necessary fraction of gravitational energy of collapse into 
kinetic energy and light of the explosive outflow. 

The basic story line of core collapse goes as follows 
(see [21 Hj for detailed reviews): At the end of a massive 
star's (8 - 10 Mq < M < 130 Mq at zero-age main se- 
quence [ZAMS]) life, nuclear burning has ceased in its 
core, which is then composed primarily of iron-group nu- 
clei (or 0-Ne nuclei at the lower end of the mass range) 
and supported against gravity by the pressure of rela- 
tivistically degenerate electrons. Eventually, the core ex- 
ceeds its effective Chandrasekhar mass and dynamical 
collapse sets in. The collapsing core separates into sub- 
sonically infalling homologous [v cx r) inner core and su- 
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personically collapsing outer core [21 [S] . When the inner 
core reaches nuclear density, the repulsive core of the nu- 
clear force leads to a stiffening of the nuclear equation of 
state (EOS). The inner core, suddenly supported by the 
stiff supernuclear EOS, overshoots its new equilibrium, 
then bounces back into the still infalling outer core. A 
shock wave forms at the sonic point between inner and 
outer core at an enclosed baryonic mass of ~O.5M0. It 
quickly moves out in radius and mass, but must do work 
in breaking up infalling iron-group nuclei. This and neu- 
trino losses from electron capture in the region behind 
the shock sap its might. The shock succumbs to the ram 
pressure of the outer core, stalls, and turns into an ac- 
cretion shock. 

The shock must be re-energized to drive a core-collapse 
supernova explosion and, in the canonical scenario, leave 
behind a neutrino-cooling and contracting protoneutron 
star. This shock revival must, depending on progenitor 
star structure, occur within ^0.5 — 3 s, otherwise accre- 
tion will push the protoneutron star over its maximum 
mass, leading to collapse and black hole formation [3]. 
Understanding the supernova mechanism, which must ro- 
bustly revive the stalled shock in supernovae from mas- 
sive stars that are observed on a daily basis, is the prin- 
ciple current challenge of core-collapse supernova theory. 

Observational clues for the supernova mechanism are 
few. Electromagnetic waves are emitted in optically thin 
regions far from the core and thus yield only second- 
hand information about the supernova mechanism. Yet, 



2 



observations of ejecta morphology, spatial distributions 
of nucleosynthetic yields and pulsar kicks are indicative 
of aspherical (i.e., multi-dimensional) processes bearing 
relevance in the explosion dynamics (e.g., [71 [S] and refer- 
ences therein). Neutrinos, on the other hand, are emitted 
deep inside the core and can provide crucial thermody- 
namic, structural, and, to some extent, dynamical infor- 
mation on what occurs in the core ^1 El 112 ■ The few 
neutrinos captured from supernova 1987A [31 [TT] have 
impressively confirmed the very basic picture of core col- 
lapse outlined in the above. 

Gravitational waves (GWs), like neutrinos emitted 
from dense regions impenetrable by photons, carry dy- 
namical information about their source. Since their emis- 
sion occurs at lowest order by accelerated quadrupole 
motions, GWs are direct probes of multi-dimensional dy- 
namics in the core that may play a crucial role in the 
supernova mechanism [3l I12j . 

Stellar collapse has long been considered a promis- 
ing source of GWs for detectors on Earth (see the his- 
torical overview in [12^ ) and much effort has gone into 
understanding the GW signature of stellar collapse and 
the subsequent evolution towards a core-collapse super- 
nova explosion. This has led to the identification of a 
range of emission processes, including rotating collapse 
and core bounce, nonaxisymmetric rotational instabil- 
ities, aspherical outflows, convection/turbulence in the 
protoneutron star and in the region immediately behind 
the shock, instabilities of the standing accretion shock, 
pulsations of the protoneutron star, asymmetric emission 
of neutrinos, and magnetic stresses (see [T31[T3] for recent 
reviews). The most recent set of simulations [12 HH- 
[19] suggest that GWs from the average core-collapse su- 
pernova may be visible throughout the Milky Way for 
the second generation of laser interferometer GW ob- 
servatories, including Advanced LIGO, Advanced Virgo, 
and LCGT [2Ql [21]. Extreme emission scenarios may 
allow detection throughout the local group of galaxies, 
including the Andromeda galaxy [T^l f33H21] , but third- 
generation detectors such as the Einstein Telescope [25] 
might hope to observe more distant events. 

If it was possible to associate an explosion mecha- 
nism with particular multi-dimensional dynamics that 
leads to a characteristic GW signal, then the detection 
or non-detection of such a signal from the next galac- 
tic core-collapse supernova could confirm or rule out this 
mechanism. To realize such an GW-observational test of 
the explosion mechanism, one must separate the signal 
from detector noise and determine its parameters, e.g., by 
matching, in some way, to signal predictions from sim- 
ulations. The most straightforward method for signal 
extraction and parameter estimation is matched filter- 
ing (e.g., ^64), which looks for a match of detector data 
with waveform templates from simulations. Matched fil- 
tering requires exact knowledge of the expected signal, 
which is possible, e.g., for the inspiral phase of compact 
binaries, since the parameter space of binary systems is 
limited, all relevant physics is understood (at least in 



the black-hole ~ black-hole binary case) and numerical 
relativity simulations can predict essentially exact wave- 
form templates (i.e., limited only by numerical error). 
However, building waveform catalogs with exact predic- 
tions for matched filtering is impossible for GWs from 
core-collapse supernovae. On the one hand, there is un- 
known physics (e.g., the nuclear EOS) and many uncon- 
strained parameters (e.g., the details of the precollapse 
configuration are poorly known) in the stellar collapse 
problem. This alone would require extensive parame- 
ter studies to build up template databases covering the 
poorly constrained parameter space. On the other hand, 
all expected GW emission processes in core-collapse su- 
pernovae are influenced or dominated by turbulent flow. 
Hence, their GW signals have a stochastic component 
that is impossible to predict, even if all initial conditions 
and physics were known exactly. Matched filtering is not 
applicable to such GW bursts. To extract the GW signal 
from the next galactic core collapse event and determine 
source physics such as the explosion mechanism, an ap- 
proach to signal extraction/reconstruction, model selec- 
tion, and parameter estimation is needed that can handle 
the stochastic nature of the expected GW signals. 

The reconstruction of both polarizations of a GW sig- 
nal requires coincident observations of two detectors; lin- 
early polarized signals can be reconstructed from data 
of just one detector and adding a third detector overde- 
termines the problem, permitting the source position on 
the sky to be determined \2!J^ . GW signal reconstruction 
was pioneered by Giirsel & Tinto (28j with a maximum 
likelihood approach, variants [551 |3D] of which have been 
incorporated into search pipelines for GW bursts [ST] . 

Summerscales et al. [32] were the first to study sig- 
nal reconstruction and parameter estimation for the GW 
burst from rotating core collapse and bounce based on 
waveforms of Ott et al. J33j. They injected signals into 
real detector noise of early LIGO science runs and used 
a maximum entropy approach to reconstruct the signal 
using data from two detectors without any a priori knowl- 
edge of the signal shape. Gross-correlation of the recon- 
structed signal with signal predictions of [33] was then 
used to determine source parameters. 

Incorporating GW signal information from core col- 
lapse simulations into detection and signal reconstruc- 
tion approaches was first considered by Brady & Ray- 
Majumder [31], who realized that the GW burst from 
rotating collapse and bounce, while being unpredictable 
in detail, has robust features that can be isolated math- 
ematically. They created an orthonormal vectorspace of 
waveforms from |331 135) using Gram-Schmidt orthonor- 
malization and isolated a subspace of essential features 
most common to all waveforms. Heng '36] also consid- 
ered waveforms from rotating core collapse and bounce 
and utilized the more recent waveform catalog of Dim- 
melmeier et al. |17j . He performed Principal Component 
Analysis (PGA; e.g., ^7\), which transforms a correlated, 
multi-dimensional data set into a set of orthogonal com- 
ponents by determining the eigenvectors and eigenvalues 



3 



of the covariance matrix of the data set. The principal 
component (PC) vectors are the eigenvectors ranked ac- 
cording to their corresponding eigenvalue, with the first 
PC being the eigenvector with the largest eigenvalue. 
Cannon et al. have also utilised this method for GW 
signals from compact binary coalescence. 

Rover et al. [35] went a step further and combined PCA 
with Bayesian inference using the Markov-Chain Monte 
Carlo technique for computing marginalization integrals 
(see, e.g., HO] for a pedagogical introduction to Bayesian 
methods) . They considered linearly-polarized waveforms 
from rotating core collapse and bounce and were able to 
reconstruct signals from modeled noise in a single detec- 
tor and infer key parameters, e.g., the nuclear EOS used 
in the simulation that led to a given trial waveform. 

In this paper, we present a proof-of-principle study 
to demonstrate that the core-collapse supernova explo- 
sion mechanism can be inferred from the GW signal of 
a galactic core-collapse supernova observed with second- 
generation GW observatories such as Advanced LIGO, 
Advanced Virgo, and LGGT. We consider the neutrino 
mechanism (e.g., [3]), the magnetorotational mechanism 
(e.g., mi), and the acoustic mechanism US US], discuss 
their essentials, and argue that they bear distinct GW 
signatures as first pointed out by Ott [121 HI] , which is 
a prerequisite for our study. We follow the approach of 
Rover et al. [39] and, for simplicity, restrict ourselves to 
a single detector, linearly polarized signals, and a Gaus- 
sian noise model at the noise level of Advanced LIGO in 
broadband mode. Like Rover et al, we adopt a Bayesian 
approach and use PCA, but, for the first time, apply it to 
multiple waveform catalogs. We associate each waveform 
catalog with one of the three mechanisms and calculate 
Bayes evidence ratios using the nested sampling algo- 
rithm [40l|45] to determine what mechanism's PCs match 
best with a given injected signal. We demonstrate that 
this approach can identify any of the considered explosion 
mechanisms with high confidence for core collapse events 
occurring at distances of up to ~2 kpc. The magnetorota- 
tional explosion mechanism can even be inferred through- 
out the Milky Way {D > 10 kpc). In addition to studying 
the explosion mechanism, we also consider the problem of 
determining whether a core collapse event is a rapidly ro- 
tating ordinary iron core collapse or an accretion-induced 
collapse (AIC) of a massive white dwarf. These two pro- 
cesses are governed by the same physics and differences 
in their waveforms are subtle [46], but our approach is 
still capable of telling them apart. 

This article is structured as follows: Section ITTl reviews 
the set of considered candidate core-collapse supernova 
mechanisms, their individual GW signatures, and the 
GW signal catalogs that we draw model waveforms from. 



In Section [nl] we introduce our method for model selec- 
tion via PCA and nested sampling. The results of our 



study are presented and discussed in detail in Section [TV] 
We summarize and conclude in Section fVl 



II. SUPERNOVA MECHANISMS AND THEIR 
GRAVITATIONAL WAVE SIGNATURES 

In this study, we consider the neutrino mechanism, the 
magnetorotational mechanism, and the acoustic mecha- 
nism for core-collapse supernova explosions and describe 
them and their characteristic GW signal features in the 
following sections. Variations of these mechanisms and 
alternatives have been discussed elsewhere (e.g., [3l I47j ) . 



A. Neutrino Mechanism 

The gravitational collapse of the iron core and the sub- 
sequent evolution of the nascent hot puffed-up protoneu- 
tron star to a cold compact neutron star release of order 
300 B (1 Bethe = 10^^ erg) of energy, ~99% of which is 
emitted in the form of neutrinos of all flavors [2] . If only 
a small fraction of the energy released in neutrinos is re- 
absorbed behind the stalled shock, leading to net heating, 
an explosion could be launched and endowed with the en- 
ergy to account for the observed range of asymptotic ex- 
plosion energies of 0.1-1 B of garden-variety core-collapse 
supernovae [48! . This is the gist of the neutrino mecha- 
nism of core-collapse supernovae, which, in its early form 
was proposed by Arnett [IS] and Colgate & White [50] . 
and in its modern form by Bethe & Wilson pTj. 

Despite its appealing simplicity, the neutrino mecha- 
nism, in its purest, spherically-symmetric (ID) form, fails 
to revive the shock in all but the lowest-mass massive 
stars with 0-Ne cores [52r55]. There is now strong evi- 
dence from axisymmetric (2D) [HirTOllSSHSD] and first 3D 
[131 [T51 simulations that the breaking of spherical 
symmetry is key to the success of the neutrino mecha- 
nism. In 2D and 3D, neutrino-driven convection in the 
region of net heating behind the shock, and the standing- 
accretion-shock instability (SASI) [57H7D] increase the ef- 
ficiency of the neutrino mechanism [60l |62l |63l [65] . 

Leaving rapid rotation aside for a moment, the 
dominant multi-dimensional GW-emitting dynamics in 
neutrino-driven core-collapse supernovae are convection 
in the protoneutron star (e.g., [75] [75]) and SASI- 
modulated convection in the region behind the stalled 
shock. GW emission from convection and SASI has been 
extensively studied in simulations in 2D [551 (TTHT^ and 
to some extent in 3D [TSl (THl jTS] [75] . The top panel in 
Fig. [l] shows a typical example waveform drawn from the 
catalog of Murphy et al. [71]. Right after core bounce, 
an initial burst of GWs is emitted by strong, so-called 
prompt convection [12] . driven by the negative entropy 
gradient left behind by the stalling shock. Subsequently, 
the GW signal settles at lower amplitudes, then picks 
up again as the SASI reaches its non-linear phase and 
high-velocity accretion downstreams penetrate deep into 
the region behind the shock, where they are decelerated, 
leading to pronounced spikes in the wave signal |71| . The 
secular rise in the signal amplitude towards the end of the 
waveform is due to the onset of an aspherical explosion 
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FIG. 1: Linearly polarized GW signal predictions for a core 
collapse event located at 10 kpc from matter dynamics in ax- 
isymmetric simulations that can be associated with the neu- 
trino mechanism (top panel, taken from [71] )i the magnetoro- 
tational mechanism (center panel, taken from [17], and the 
acoustic mechanism (bottom panel, taken from [T^). Note 
the varying ranges of the time and strain axes. See text for 
discussion. 



[T5|, 1581 [71] , but occurs at too low characteristic frequen- 
cies to be visible to Advanced LIGO-class detectors. Not 
included in the top panel of Fig. [T] is the contribution to 
the GW signal from anisotropic neutrino emission [771 - 
179] ■ which can dominate in amplitude, but, like the con- 
tribution from aspherical outflow, occurs on timescales 
too long to lead to emission at frequencies detectable 
by the upcoming generation of ground-based detectors 
[l5l[l6l[Ml[72l[74l[76]. 

Overall, the detectable GW signal from a neutrino- 
driven nonrotating or slowly rotating core-collapse super- 
nova will have random polarization, a broadband spec- 
trum from ~100 — 1000 Hz and typical strain amplitudes 
\h\ of order 10~^^ (D/10 kpc)~^, with individual peaks 
reaching lO'^i (D/IO kpc)"! [H [H]. The typical du- 
ration of emission is the time from core bounce to the 



launch of the explosion, 0.3 — Is, but convection inside 
the cooling protoneutron star can continue to emit GWs 
at lower amplitudes and higher frequencies for many sec- 
onds afterwards [TH [75] . Typical total emitted GW en- 
ergies are in the range 10"" - 10"^ Mq(? [ni[551[7T] . 

The effects of rotation on the neutrino mechanism and 
its GW signature are not yet fully understood (see, e.g., 
[Tl[ini[ll|S71[Ml[Eni[HI]) and it cannot be excluded that 
contributions from rotational dynamics may modify the 
GW signal of neutrino-driven core-collapse supernovae. 
However, results from the systematic rotating core col- 
lapse studies of [TT] [ISl HIl HI] suggest that once rotation 
rates become sufficiently high to alter the dynamics, the 
explosion is actually more likely to occur via the magne- 
torotational mechanism discussed in §H B[ This, however, 
is under the provision that the magnetorotational insta- 
bility (e.g., [HH IH3]) works robustly and builds up the 
required strong magnetic fields to drive an explosion. 

Keeping the above caveats in mind, for the purpose of 
this study, we make the assumption that the GW signa- 
ture of neutrino-driven core-collapse supernovae is unaf- 
fected by rotational effects. 



1. GW Signal Catalogs 

In this study, we use the catalog of Murphy et al. [7T] , 
which is available for download from [84]. The Mur- 
phy et al. catalog (in the following, we refer to wave- 
forms from this catalog as Mur waveforms) encompasses 
16 waveforms that were extracted via the quadrupole 
formula (e.g., [85]) from Newtonian axisymmetric core 
collapse simulations that used a parameterized scheme 
for electron capture and neutrino heating/cooling and 
included only the monopole component of the gravita- 
tion potential as described in [60] [71] . The Murphy et al. 
simulations are nonrotating and the parameter space cov- 
ered is spanned by progenitor ZAMS mass ({12, 15, 20, 
and 40} Mq ) and by the dialed- in total electron and anti- 
electron neutrino luminosity. 

Yakunin et al. [58] performed self-consistent axisym- 
metric Newtonian (with an approximate-GR monopole 
term of the gravitational potential |86| ) radiation- 
hydrodynamics simulations of neutrino-driven core- 
collapse supernovae. They provide three waveforms at 
[57] , obtained from simulations using progenitors of (12, 
15, and 25) Mq. We use the Yakunin waveforms (de- 
noted, in the following, as Yak waveforms) to test the 
robustness of our supernova mechanism determination 
algorithm, which uses the PCs of the Mur waveforms. 

Since we are limiting ourselves to one detector in this 
proof-of-principle study, we are considering only linearly 
polarized signals. Gravitational waveforms with -I- and 
X polarizations from 3D simulations of neutrino-driven 
core-collapse supernovae [151 116[ 176) will be considered in 
future work. 
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B. Magnetorotational Mechanism 

The conservation of angular momentum in core col- 
lapse to a protoneutron star leads to a spin-up by a fac- 
tor of ~1000 [HB]. Starting from a precollapse angular 
velocity distribution that may be expected to be more 
or less uniform in the inner core (e.g., [89]), homologous 
collapse preserves the uniform rotation of the inner core 
while the supersonic collapse of the outer core leads to 
strong differential rotation in the outer protoneutron star 
and in the region between protoneutron star and shock 

A rapidly spinning precollapse core with a period of 
order 1 s results in a ms-period protoneutron star, with a 
rotational energy of order 10 B, which is about ten times 
greater than the typical core-collapse supernova explo- 
sion energy. If only a fraction of this energy was tapped, 
a strong explosion could be triggered. 

Theory and simulations (e.g., [H EH lOHS]) have 
shown that magnetorotational processes are efficient at 
extracting spin energy and can drive collimated outflows, 
leading to energetic bipolar jet-like explosions. Recent 
work [IHl suggests that magnetic fields of the 

order of 10^^ G with strong toroidal components are re- 
quired to yield the necessary magnetic stresses to drive 
a strong bipolar explosion. If 10^^ G fields were to arise 
from flux compression in collapse alone, precollapse core 
fields would have to be of order 10^^ G [HI which 
is about 3 orders of magnitude larger than predicted by 
stellar evolution models (e.g., [521117]). It is more likely 
that the most significant amplification occurs after core 
bounce via rotational winding of poloidal into toroidal 
field (a linear process) , the non-linear magnetorotational 
instability (MRI, which is not yet fully understood in the 
core collapse context [5^ 183]). Both processes operate 
on the free energy stored in differential rotation, which 
is abundant in the outer core. 

For the magnetorotational mechanism to work, pre- 
collapse spin periods < 4 — 5s appear to be required 
[4T] . Such rapid rotation leads to a strongly centrifugally- 
deformed inner core with a large quadrupole moment 
{£ = 2; due to its oblateness), which rapidly changes 
during core bounce, leading to a strong burst of GWs. 
The GW signal from rotating collapse and bounce has 
been studied extensively and the most recent general- 
relativistic simulations have shown it to be of rather 
generic morphology with a single strong peak at bounce 
and a subsequent ringdown as the protoneutron star core 
settles into its new equilibrium [TT] [HHl [H]. A typi- 
cal example GW signal taken from the catalog of Dim- 
melmeier et al. [TT] 1100] is shown in the center panel of 
Fig. [Tj The core collapse and bounce phase proceeds es- 
sentially axisymmetrically even in very rapidly spinning 
cores [m 11011 1102] and its GW signal is linearly polar- 
ized with vanishing amplitude seen by an observer lo- 
cated along the symmetry axis and maximum amplitude 
for an equatorial observer. Typical emission durations for 
the linearly polarized GWs from core bounce are of or- 



der 10 ms and peak GW amplitudes for rapidly spinning 
cores that may lead to magnetorotational explosions are 
of order 10~^^ — 10^^" at lOkpc with most of the energy 
being emitted around 500 — 800 Hz in cores that reach 
nuclear density and bounce due to the stiffening of the 
nuclear EOS. Cores with initial spin periods shorter than 
^0.5 —Is experience a slow bounce at sub-nuclear densi- 
ties strongly influenced or dominated by the centrifugal 
force. They emit most of the GW energy at frequencies 
below ~ 200Hz [ni in]. Typical emitted GW energies 
are in the range 10-^° - IQ-^MqC^. The GW signal 
from rotating collapse and core bounce is unlikely to be 
affected by MHD effects, since the build up to dynam- 
ically relevant field strengths occurs only after bounce 

[iHiiiiiiniissj- 

Due to the strong rotational deformation of the pro- 
toneutron star, neutrinos decouple from the matter at 
smaller radii and hotter temperatures in polar regions 
than near the equator. This leads to the emission of 
a larger neutrino flux with a harder neutrino spectrum 
in polar regions (e.g., [59] )• This globally asymmet- 
ric neutrino emission results in a secularly rising low- 
frequency GW signal |12] . Similar low- frequency contri- 
butions will come from the bipolar outflow characteris- 
tic for a magnetorotational explosion and from magnetic 
stresses [13 |M) . The low- frequency waveform com- 
ponents are not shown in the center panel of Fig. [l] and 
are not detectable by the upcoming second-generation 
earthbound GW observatories. 

Also associated with rapid rotation and the magnetoro- 
tational mechanism are rotational instabilities that may 
lead to nonaxisymmetric deformations of the protoneu- 
tron star whose "bar-mode" (to — 2) components may 
emit elliptically polarized GWs for tens to hundreds of 
milliseconds [Ml I101H104] . However, these instabilities, 
and in particular their interplay with magnetic fields and 
the MRI (see, e.g., [105] ). are not yet fully understood. 
Since we are considering only linearly polarized signals 
and are limited to one detector, we do not include GW 
signals from these nonaxisymmetric instabilities in this 
study. 



1. GW Signal Catalogs 

We employ the large (128 waveforms) GW signal cat- 
alog of Dimmelmeier et al. [TT] 1100] (Dim in the follow- 
ing) , who performed 2D GR simulations of rotating iron 
core collapse for (11.2, 15, 20, and 40) Mq progenitors 
and two different nuclear EOS, varying initial rotation 
rate and degree of differential rotation. They approxi- 
mated the effects of electron capture during collapse by 
parametrizing the electron fraction Yg as a function of 
density, which yields inner core sizes that are very close 
to those obtained with full neutrino transport [106] . The 
inner core size determines the amount of mass and angu- 
lar momentum that can be dynamically relevant during 
core bounce and, hence, is a determining factor in the 
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GW signal [55] • The Dim catalog was also used by the 
previous parameter estimation work of Rover et al. |39j . 
For testing, we use the three additional Dim waveforms 
computed for 139] that are not part of the original Dim 
catalog and were used to test their algorithm. We label 
this set of extra waveforms as DimExtra. 

For studying the robustness of our mechanism- 
determination approach, we draw gravitational wave- 
forms of rotating models from the catalog of Scheideg- 
ger et al. [10211107] . (Sch in the following) who performed 
3D Newtonian-MHD rotating iron core collapse calcula- 
tions with a spherical approximate-GR gravitational po- 
tential and employed the same EOS and electron capture 
treatment as Dimmelmeier et al. [17j . but used different 
progenitor models. 

Furthermore, we use the GW signal catalog of Abdika- 
malov et al. [13 [51] (Abd in the following) who used the 
same numerical code as Dimmelmeier et al. [17' , but stud- 
ied the rapidly spinning accretion-induced collapse (AIC) 
of massive white dwarfs to neutron stars. This process 
yields a GW signal very similar to rotating iron core col- 
lapse and explosions in AIC may occur also via the mag- 
nctorotational mechanisms [108j . We include this catalog 
of 106 waveforms to see if our algorithm can differenti- 
ate between rotating iron core collapse and rotating AIC 
assuming the Dim and Abd catalogs correctly predict the 
respective GW signals. 



C. Acoustic Mechanism 

The core-collapse supernova evolution in the acoustic 
mechanism proposed by Burrows et al. [551 IHl SS] is ini- 
tially identical to the one expected for the neutrino mech- 
anism. Neutrino heating, convection and the SASI set 
the stage, but no explosion is triggered for >500 ms after 
bounce. At this point, the SASI is in its highly non- linear 
phase and modulates high-velocity accretion downflows 
that impact on the protoneutron star and excite core pul- 
sations (primarily £ = {1,2} g-modes). Over hundreds of 
milliseconds, these pulsations reach large amplitudes and 
damp via the emission of strong sound waves. Traveling 
down the steep density gradient in the region behind the 
shock, the sound waves steepen to shocks and dissipate 
their energy behind and in the shock. This mechanism is 
robust in the simulations by Burrows et al. [52 Ull [13] , 
but requires >ls to develop, thus leads to massive NSs, 
and tends to yield explosion energies on the lower side of 
what is observed. 

The GW signature of the acoustic mechanism is dom- 
inated by the strong emission from the quadrupole com- 
ponents of the protoneutron star core pulsations that are 
quasi-periodic (their frequency shifts secularly along with 
the changing protoneutron star structure) and become 
very strong >800 — 1000ms after core bounce [121 122] . 
The lower panel of Fig. [T] depicts a typical example wave- 
form from Ott et al. [Il[21], who studied the GW sig- 
nature of the acoustic mechanism based on the simula- 



tions of Burrows et al. [121 US]. At early times, the GW 
signal is essentially the same as expected for the neu- 
trino mechanism, but once the protoneutron star core 
pulsations grow strong, they are hard to miss. The sim- 
ulations of Burrows et al. [35] HH US] were axisymmet- 
ric and the resulting GW signals are linearly polarized, 
though in 3D, one would expect oscillation power also in 
non-axisymmetric components. Typical maximum strain 
amplitudes are of order few x 10~^^ — 10~^° and mul- 
tiple modes with frequencies between ~ 600 — 1000 Hz 
contribute to the emission. Since the pulsations last for 
many cycles, the emitted GW energies may be large and 
are predicted to be of order 10~* — 1O~^M0C^ and ex- 
treme models reach few x IO^^A/qC^ [T2ll22]. 

There are multiple caveats associated with the acoustic 
mechanism that must be mentioned. Most importantly, 
the acoustic mechanism has been found in simulations of 
only one group with a single simulation code, but oth- 
ers have not yet ruled out the possibility of strong pro- 
toneutron star pulsations at late times (e.g., [19j). In 
a non-linear perturbation study, Weinberg & Quataert 
[109] found that the protoneutron star pulsation ampli- 
tudes may be limited by a parametric instability involv- 
ing high-order modes that damp efficiently via neutrino 
emission and are not presently resolved in numerical sim- 
ulations. This would limit the protoneutron star pulsa- 
tions to dynamically insignificant amplitudes. Moreover, 
the simulations of Burrows et al. were axisymmetric and 
nonrotating or only very slowly rotating. It is not clear to 
what amplitudes individual protoneutron star pulsation 
modes would grow in 3D. Rapid rotation, due to its sta- 
bilizing effect on convection and SASI [44] [59] , may likely 
inhibit the growth of pulsations. Both 3D and rotational 
effects remain to be explored. 



1. GW Signal Catalogs 



We employ the set of 7 waveforms from the models of 
[13] analyzed by Ott [T5] and available at [51]. We refer 
to this set as the Ott catalog in the following and use 
them to compute PCs for the acoustic mechanism's GW 
signature. All waveforms were computed on the basis of 
the Burrows et al. [55] US] US] simulations and differ only 
in the employed progenitor model, covering a range in 
ZAMS mass from 11.2 to 25 Mq. 

Three additional waveforms of an earlier study of 
Ott et al. [22] are available [84]. We label this small set 
OttExtra and use them for testing our method's capa- 
bility of correctly identifying them as coming from stars 
exploding via the acoustic mechanism. 
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III. DATA ANALYSIS 



A. Strategy 



The three example gravitational waveforms shown in 
Fig.[T]that are associated with the three supernova mech- 
anisms are clearly different. Provided the assumptions 
made in associating these signals with the various mech- 
anisms are correct, a GW signal detected from a core- 
collapse supernova should, in principle, allow to deter- 
mine the explosion mechanism. To do so in practice, two 
problems most be overcome: (i) The exact waveform of 
an incident signal is impossible to predict in advance. 
(ii) Real GW detectors are noisy instruments (see, e.g., 
[110) for a discussion of detectors and noise sources) and 
any GW signal will be contaminated by detector noise. 
In other words, it is necessary to develop a data anal- 
ysis algorithm that is capable of distinguishing between 
underlying physical models (e.g., supernova mechanisms) 
on the basis of a noisy signal whose detailed shape cannot 
be predicted exactly. 

In the following subsections, we describe the compo- 
nents of a Bayesian data analysis algorithm which classi- 
fies detected GW signals from core-collapse supernovae as 
belonging to one of a set of signal catalogs, representing, 
e.g., different explosion mechanisms. A block diagram 
of the analysis algorithm, which we call the Supernova 
Model Evidence Extractor (SMEE), is shown in Fig. [2] 
SMEE is implemented in MATLAB^ 

In a first step, SMEE performs Principal Component 
Analysis (PCA) via singular value decomposition (SVD) 
on the waveforms in each catalog to create sets of or- 
thogonal basis vectors, the Principal Components (PCs), 
which are ordered according to their prevalence in their 
catalog. In other words, the first PC represents the most 
common feature of all signals in the catalog, the second 
PC represents the second most common feature, and so 
on. Using a complete set of PCs, each waveform can be 
reconstructed as a linear combination of PCs for the cor- 
responding catalog, allowing each waveform to be simply 
parameterized by the PC coefficients in the linear com- 
bination. However, since PCs are expected to span the 
parameter space defined by each catalog of waveforms 
efficiently, catalog waveforms may be reconstructed with 
good accuracy already with a set of PCs that is signifi- 
cantly smaller than the number of waveforms in the cata- 
log. Moreover, non-catalog waveforms (i.e., real signals) 
may be identified as belonging to the same class of sig- 
nals as catalog waveforms if they can be approximately 
matched with the first few PCs of a catalog. 

SMEE then uses Bayesian model selection and com- 
putes the logarithm of the Bayes factor to distinguish be- 
tween GW signal classes. The Bayes factor is the ratio of 
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FIG. 2: Block diagram of the Supernova Model Evidence Ex- 
tractor (SMEE). A desired core-collapse supernova gravita- 
tional wave signal is injected into noise, and the algorithm 
compares it to the principal components (PCs) of a given 
waveform catalog representing a particular model. The PCs 
are constructed via singular value decomposition (SVD). The 
sign of the log Bayes factor between two PC sets indicates 
which model is favored by the data. 



the evidences for two competing hypotheses and, for the 
purpose of our analysis, we weigh the evidence that the 
observed data supports the presence of a GW signal con- 
sistent with signals from one of any two competing cat- 
alogs. This requires summing up the likelihood function 
times the prior across all possible signal parameters (in 
our case, values of PC coefficients) to determine the evi- 
dence (also called the marginal likelihood) for two differ- 
ent signal models to be tested. SMEE accomplishes this 
efficiently via the Nested Sampling algorithm ^40, ^45^ . 



B. Bayesian Model Selection 



1 The MathWorks Inc.. 



Natick, 



MA 01760, 



Ihttp : //www . mathworks . com/products/matlab/ 



USA. In our analysis, we employ Bayesian Model Selection, 

similar to that described in [lllj . Specifically, we use the 
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Bayes factor to compare the probabilities of two compet- 
ing models. In general terms, the Bayes factor Bij can 
be written as the ratio of the evidences p{D\M), 



B, 



p{D\M,) 



(1) 



where Mi and Mj are two competing models tested us- 
ing the data D. The evidence for each model is obtained 
by integrating the product of its likelihood function and 
prior across all model parameter values 9, such that, 



p{D\M)= p{e\M)p{D\e,M)d9 . 



(2) 



The evidence will be greater for a model that is sup- 
ported by the data. Therefore, the Bayes factor indicates 
which of the two competing models is preferred by the 
data. It is often more convenient to compare models by 
using the natural logarithm of the Bayes factor. 



log By - logp{D\M,) - \ogp{D\Mj 



(3) 



In this case, log By > means Mi is the preferred 
model whereas log Bij < will point to Mj being favored. 



distribution is only found as a by-product. To find the 
evidence, a set of "live points" are found through creating 
a stochastic sampling of the prior distribution to gener- 
ate a set of N samples which are denoted as di, where 
i = 1 ... TV. The evidence integral (Eq. [2]) is then written 
as 



p{D\M)^ p{e\M)p{D\e,M)de , 



N 



1=1 

N 

^ LiWi , 



1=1 



where the weight. 



(4) 



(5) 



is the fraction of the prior distribution represented by 
the z-th sample and Li is its likelihood. It is this weighted 
likelihood that is calculated by the Nested Sampling al- 
gorithm and subsequently used to obtained the evidence. 
More details on the Nested Sampling algorithm can be 
found in [451ITT2]. 



C. Nested Sampling Algorithm 

For evaluating log Bij, wc first need to calculate the 
evidences p{D\Mi) and p{D\Mj) for the two models Mi 
and Mj. From Eq. [2] we see that the evidence is the sum 
of the likelihood times the prior determined for all pos- 
sible parameter values of the desired model. An exhaus- 
tive, brute-force approach to computing the evidence by 
calculating the likelihood values for every choice of pa- 
rameter values is computationally prohibitive. It is also 
an inefficient way of determining the evidence since the 
likelihood values will be most significant, and therefore 
contribute most to the evidence, for a small subset of pa- 
rameter values which constructs a waveform that closely 
resembles the data. For most other combinations of the 
model's parameters, the likelihood will be insignificant 
and not contribute to the evidence. 

Therefore, we choose to follow the approach of Veitch 
et al. |112j and employ Nested Sampling [301 US] to effi- 
ciently calculate the evidence integral. The Nested Sam- 
pling algorithm determines the evidence integrals by cal- 
culating the likelihood for a selected sample of param- 
eter values for the desired model. Initially, the model's 
parameter values are randomly selected before the algo- 
rithm iteratively converges on the set of parameter val- 
ues that produce the most significant likelihood values. 
It is similar to the Markov-Chain Monte Carlo approach 
(e.g., [40]) except that the primary output of the Nested 
Sampling algorithm is the evidence, which can be im- 
mediately obtained by summation whereas the posterior 



D. Principal Component Analysis via Singular 
Value Decomposition 

Each core-collapse supernova waveform catalog con- 
sists of a number of GW signals obtained for different 
initial conditions and simulation parameters (e.g., pro- 
genitor star mass, EOS, rotational configuration etc.). 
While individual waveforms of one catalog are different 
in detail, they generally exhibit strong common general 
features. This can be exploited by principal component 
analysis (PCA) [37], which isolates the most common 
features of waveforms in linearly independent principal 
components (PCs) ordered by their relevance. The first 
few PCs may already be sufficient to efficiently span their 
entire catalog, as was shown in [3S] and [35] for the Dim 
catalog (see Sect. |IIB"T| ). 

The PCs are obtained via singular value decomposition 
(SVD) (e.g., [37]) of time-domain waveforms from each 
catalog. To perform SVD on a catalog with m waveforms, 
a matrix A is created such that each of its columns cor- 
responds to a waveform of uniform length n from the 
catalog. 

The n X m waveform matrix A is factored so that 

A = USV^, (6) 

where U is an n x n matrix whose columns correspond 
to the eigenvectors of AA^ . Similarly, the columns of 
the m X m matrix V correspond to the eigenvectors of 
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A'^A and S is an n x m diagonal matrix whose elements 
correspond to the square root of the corresponding eigen- 
values. 

Since AA"^ is the covariance matrix of A, the eigen- 
vectors in U are effectively an orthonormal basis which 
span the m-dimensional parameter space defined by the 
catalog of waveforms used to construct A. Note that, in 
practice, n ^ m and it is impractical to determine the 
eigenvectors in U directly. Instead, the smaller V and 
its corresponding eigenvalues in S are first determined 
which are subsequently used to derive U. 

The orthonormal eigenvectors of U are the PCs and 
are ranked by their corresponding eigenvalues. The PC 
with the largest corresponding eigenvalue is referred to 
as the first PC and consists of the most significant com- 
mon features of all waveforms in the catalog. It follows 
that the PC with the second largest corresponding eigen- 
value is the second PC and consists of the second most 
significant common features and so on. 

The waveforms in A can be reconstructed by taking a 
linear combination of PCs, 



(7) 



integration in Eq. [2] and obtain an analytic form for the 
noise evidence function, 



N 



p{D\M„)^T\ =exp 



i=l 



(9) 



In both Eq.[8]and|9] the standard deviation of the noise 
is a function of each sample in the data, because the sim- 
ulated noise is designed to correspond to the expected 
sensitivity of Advanced LIGO, which varies as a function 
of frequency (see Sec. IIIF). To handle the frequency- 



dependent noise, the signal and noise evidences are calcu- 
lated in the frequency domain with Di, hi(/3) and tXi cor- 
responding to data, reconstructed waveform, and noise 
in the i-th frequency bin respectively. In particular, each 
reconstructed waveform is obtained by taking a linear 
combination of the Fourier transforms of its correspond- 
ing PCs. 

The natural logarithm of the Bayes factor used to com- 
pare the signal model to the noise model is then simply 



log BsN = log[p(L'|M,)] - log[p{D\M^)]. 



(10) 



where hi is the desired waveform from the catalog, Uj 
is the jth PC from the U matrix and f3j is the corre- 
sponding PC coefficient, which can be obtained by pro- 
jection of hi onto Uj. The sum of k PCs produces an 
approximation of the desired waveform since k < m. 



E. Signal and Noise Models 

For the analysis described here, two types of models 
are considered. The signal model Ms tests the presence 
of a signal waveform h{f3) in the data. Here, PCA is per- 
formed for each catalog (using SVD; see section HID) 



and each waveform is parameterized by its PC coefh- 
cients (/?). The Gaussian likelihood function for the sig- 
nal model is 



N 



(A 



2a? 



, (8) 



where cr^ is the standard deviation of the noise, /ii(/3) 
is the desired waveform reconstructed from the PCs and 
N is the length of the data with a corresponding index 
i. The evidence for the signal model is determined by 
performing the integral in Eq. [2] numerically across all 
chosen values of (3. 

On the other hand, the noise model Af„ tests the data's 
consistency with Gaussian noise. The likelihood function 
for the noise model is the same as that in Eq. [8j but with 
/i(/3) = 0. From this, it is straightforward to perform the 



F. Generation of Simulated Noise 

We generate Gaussian colored noise, assuming a sin- 
gle Advanced LIGO detector in the proposed broadband 
configuration (the so-called "zero detuning, high-power" 
mode). We employ the data file ZERO_DET_high_P.txt 
provided by |113j , which contains \JS{f), the square root 
of the one-sided detector noise power spectral density in 
units of (Hz)"^/^. An open-source implementation for 
MATLAB of what we describe in the following can be found 
in [TTi] . 

The real discrete time-domain noise n(tj), where tj de- 
notes the j-th discrete time interval of size Ai, is obtained 
by inverse discrete Fourier transform from the complex 
frequency-domain noise n{fk), where fk denotes the k- 
th discrete frequency interval of size A/ = l/{NtAt), 
where Nt is the number of intervals in the time domain. 
Since the time-domain noise is real and has zero mean, 
the frequency-domain noise must obey 

n{-.f) = n*if) , (11) 

h{,f = 0) = nif = /Nyq) = Hf = -/Nyq) - 0.(12) 

Here, /Nyq — l/(2Ai) is the Nyquist frequency. If Nt is 
the even number of equally spaced bins in time of width 
At, then Nf = Nt/2 — 1 is the number of independent 
frequency bins fk in the frequency domain of width A/. 
The frequency variable fk assumes values from —/Nyq to 

/Nyq- 

l^-l/fc)! = V n{fk)n*{fk) is a two-sided amplitude spec- 
tral density. We generate n{fk) by sampling the standard 
normal distribution (zero mean, variance one) weighted 
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FIG. 3: The first tliree Principal Components (PCs) of tlie waveforms from the Dimmelmeier (Dim) [17lllUUj . Murphy (Mur) |71l 
I84| . and Ott |121I84| catalogs, which we take to be representative of the magnetorotational, neutrino and acoustic mechanisms, 
respectively (see Sees. II A |II B and II C I. In creating the PCs, the waveforms of each catalog are placed in a systematic way 
in a 3 s interval and padded left and right by zeros. The vertical axis is a dimensionless scale which represents the amplitude, 
which we have normalized by the maximum amplitude over all 3 PCs of each catalog shown here. 



by the noise transfer function T(/fc) = ^JS{fk)/V2. The 
real and imaginary parts of n(/fe) are then given for each 
fk e (0, /Nyq) by 



mfk)) 



1 

1 



T(/) RANDN , 
T(/) RANDN , 



(13) 
(14) 



where RANDN is a random number sampled from the stan- 
dard normal distribution. The remaining h{fk) are then 
obtained via Eqs. [TT] and [T2j 

The inverse discrete Fourier transform oin(fk) to time 
domain noise n{tj) will preserve its Gaussian character 
in the time domain in the limit of small sampling inter- 
val [39l I112j . Specifically, we use the following definition 
of the discrete Fourier transform (DFT) for transform- 
ing noise from the frequency to the time domain when 



needed: 



7V-1 



nifk) = J2 «(ij)exp(-27rijfc/iV) 



= ^ E "(/fc)exp(2^ij/c/iV) 



(15) 



(16) 



k=0 



For convenience, we define the matched filter signal- 
to-noise ratio (SNR) of a GW signal h as 



SNR2 = 4 



Mf) 



Sif) 



-df 



2 ^ r^'^) 



(17) 
(18) 
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where 5'(/) is the one-sided noise power spectral density. 
The factor At^ is applied to correct the dimensions of 
h{f), which we obtain via the DFT defined by Eq. 



15 



G. Application in SMEE 

1. GW Signal Preparation and PC A 

Before carrying out PCA and injecting signals into 
noise, all waveforms are buffered with zeros to be of 
length n, which we choose to correspond to 3 s at a sam- 
pling rate of 4096 Hz, allowing us to comfortably accom- 
modate the longest available core-collapse supernova GW 
signals. The Advanced LIGO sampling rate is 16 kHz. 
The reduced sampling rate we choose saves computation 
time and is sufficient to capture the frequency content of 
the core-collapse supernova waveforms considered here, 
which have most of their power at ~50 — 1000 Hz. 

We align waveforms from the Dim catalog at their max- 
imum (the spike at core bounce). Waveforms from the 
Mur and Ott catalogs are aligned so that the onsets of 
emission coincide. All waveforms are shifted so that they 
are aligned to the 4000-th point in the SMEE input data 
file, corresponding to about the 1 s mark in the 3 s inter- 
val, to leave ample space left and right of the waveform. 

In Fig. [3) we present the first three PCs computed for 
the Dim (magnetorotational mechanism; left panel), Mur 
(neutrino mechanism; center panel), and Ott (acoustic 
mechanism; right panel) catalogs. Before generating PCs 
for the Mur catalog we filter out the secular low-frequency 
drifts present in the Mur waveforms (see Fig. [l]) by high- 
passing the signal above 30 Hz. Since the low-frequency 
components are hidden in detector noise even when the 
source is nearby, dropping them improves the efficiency of 
our subsequent Bayesian analysis and signal reconstruc- 
tion. We apply the same high-passing to trial waveforms 
for the neutrino mechanism before injecting them into 
noise. 



2. Signal Injection and Model Selection 

We inject trial GW signals into simulated Gaussian 
Advanced LIGO noise and use SMEE to determine which 
signal model (e.g., what core-collapse supernova explo- 
sion mechanism) a given injected signal belongs to via 
the evaluation of the logarithmic Bayes factors log Bsn 
(Eq. 10) for an injected signal for each signal model S 
and the noise model N. Comparing two signal models 
i and j is then accomplished by computing log Bij = 

logB^N - log BjN- 

SMEE's model selection operates in the frequency do- 
main. Trial GW signal and the PCs belonging to the sig- 
nal model under consideration are transformed into the 
frequency domain via DFT and the trial GW signal is 
added to the complex frequency-domain noise, retaining 
phase information. The Nested Sampling algorithm is 



then invoked to marginalize the PC coefficients (3^. The 
prior for each coefficient is flat and uniform. The prior 
range for each (3k is determined by first reprojecting all 
waveforms of a given catalog back onto the PCs to com- 
pute /3ki for each PC k and waveform I of the catalog. 
The range of expected possible values of f3k is then found 
by taking the minimum and maximum of jSki over all I 
and adjusting these numbers by 10% down and up to add 
a margin of error to account for uncertainty due to the 
noise, motivated by the findings of f55^. 

Keeping the noise model fixed, the results of SMEE's 
computations will depend on the SNR of the signal, i.e. 
the distance to the core collapse event, and on the amount 
of information we can provide to SMEE about expected 
signals in the form of PCs. 

The maximum number of PCs at SMEE's disposal is 
limited by the number of waveforms used to determine 
the set of PCs. While each catalog used in this study has 
a different number of waveforms, we choose to simplify 
our analysis by using the same number of PCs for all cat- 
alogs. Hence, the maximum number of PCs we use here is 
7 and is set by the number of waveform in the Ott catalog 
(see |IIC 1| . Using 7 PCs gives SMEE complete infor- 
mation about signals belonging to the Ott catalog and 
significant, but incomplete information about waveforms 
from the Dim, Mur, and Abd catalogs. We also carry out 
SMEE runs with less than 7 PCs to study the dependence 
on the number of PCs employed. Using only a small sub- 
set of a catalog's PCs limits SMEE's ability to precisely 
reconstruct injected catalog waveforms, but it represents 
the real-life situation that the a priori information about 
a detected signal is severely limited. Our goal here is not 
to ideally reconstruct signals but to show that determin- 
ing the underlying physical model of an observed signal 
is possible with limited advance knowledge. 



IV. RESULTS 

A. Response to Gaussian Noise 

For interpreting the results of SMEE's Bayesian model 
selection on the basis of Eq. [3) it is necessary to quan- 
tify and understand SMEE's response to pure Gaussian 
detector noise without a signal being present. To this 
end, we run SMEE on 10,000 randomized instances of 
Advanced LIGO detector noise (generated as described 
in pIIF ) without injecting signals and compute log Bsn 
(Eq. |10[) in the absence of a signal for each signal model 
S. The results, shown in Fig. |4j follow a Gaussian distri- 
bution with a mean corresponding to the expected value 



— J2iLi ^2a'^ ; where N is the number of PCs employed 
(see j |IlIE |. The average logarithmic Bayes factors ob- 
tained for 10,000 instances of noise indicate that noise, 
or any signal fully consistent with noise, is most likely to 
have a logarithmic Bayes factor of —54.0 when SMEE 
is run with 7 PCs of the Dim catalog. For the Ott, 
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FIG. 4: Results from running SMEE with 7 PCs of the Dim 
catalog and without an injected signal on 10,000 randomized 
instances of Gaussian Advanced LIGO noise, generated as 
described in Sec. |III F] A signal consistent with noise is most 
likely to have a logarithmic Bayes factor of ~— 54.0. The red 
line plots a Gaussian distribution with a mean of —53.96 and 
a standard deviation a = 0.17. 



Mur, and Abd catalogs, the expected logarithmic Bayes 
factors for pure Gaussian noise and 7 PCs are —52.1, 
—52.3, and —53.0, respectively. The observed expecta- 
tion value are very comparable to those calculated for the 
Dim (-53.9), Dtt (-52.2), Mur (-52.3), and Abd (-52.9) 
catalogs, respectively, verifying that SMEE is operating 
as expected. We have repeated this experiment for the 
case when only 3 PCs are used and also in this case find 
that SMEE closely reproduces the predicted expectation 
values, which are near —26 in the 3-PC case. 

Since the logarithmic Bayes factors follow a Gaussian 
distribution, we can set a threshold using the standard 
deviations as an indicator for the expected false alarm 
rate. Ideally, for the Dim catalog, a 1% false alarm rate 
would correspond to a threshold that is ^-^2.6 times the 
standard deviation, corresponding to ~0.44 above the 
mean. However, we note that the expected logarithmic 
Bayes factor value varies between different catalogs and, 
for a fixed false alarm rate, we would require a different 
threshold for each catalog. This variation can be ad- 
dress by re-normalising all Bayes factors so that they are 
the same for all catalogs when there is only noise. But, 
since the focus of our work here is to distinguish between 
different waveforms and not to perform a study on the 
detection efficiency of GW signals, we choose to take the 
more conservative approach of simply setting a higher 
threshold. Therefore, we conservatively choose to iden- 
tify a signal as being distinct from noise if its log Bsn is 
greater than —47 (in cases in which we use 3 PCs, this 
number is —21). When comparing two signal models Mi 
and Mj, we conservatively identify model Mi as favored 
if log Bij > 5 (and vice versa). 
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FIG. 5: Mean logBsjv as a function of signal-to-noise ratio 
(SNR; Eq. |18[ ) for 5 representative waveforms from the Mur, 
Ott and Dim catalogs using 7 Principal Components (PC). 
The shaded areas represent the standard error in the mean 
value of logBsAT for each waveform catalog computed as ct = 
±7V~^(Ei(2; — Xi)^Y^'^ , where x is the mean and Xi are the 
individual SNRs and A*" is the number of waveforms. Values 
of logSsjv below —47 in the 7-PC ase and below —21 in 
the 3-PC case indicate that the algorithm considers it more 
likely that there is no signal detectable in the noise. Table |l] 
summarizes numerical results for the minimum SNR for which 
log Bsn > -47. 



B. Signal vs. Noise 

The minimal GW signal strength required for SMEE to 
be able to select the core-collapse supernova mechanism 
is an important question. The primary prerequisite for 
an incident GW signal to be useful for model selection 
is that SMEE can distinguish it from detector noise, i.e., 
we must find the minimum signal strength (i.e., SNR) so 
that log Bsn > —47 (when 7 PCs are used; Eq. 10 and 
SjlVA|. 



In order to determine the range of minimum SNR re- 
quired across and within core-collapse supernova GW sig- 
nal types, we draw 5 representative waveforms from the 
Dim, Mur, and Ott catalogs and run them through SMEE 
at varying SNR, using 7 PCs generated from the catalog 
to which each injected waveform belongs. The result of 
this exercise is shown in Fig. [5] and summarized in Tab.jl] 
Generally, an SNR > 4 - 5 is required for SMEE to find 
log Bsn > —47 in the idealized setting that we consider 
here. In a real fully blind search, unknown arrival times 
and non-Gaussianity of real detector noise will generally 
require an SNR in excess of 8 for a detection statement 

(e.g., m)- 
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TABLE I: Minimum SNR at which the injected waveform will 
give a log Bsn above the threshold defined in §IV A| This is 
shown for 5 representative waveforms from the Dim, Mur, and 
Ott catalogs, which we take to be representative of the magne- 
torotational, neutrino, and acoustic mechanisms, respectively. 
We provide results from SMEE runs with 3 and 7 PCs. 



Waveform 
Name 


Minimum SNR 
3 PCs 7 PCs 


Dim [17] 






sllaloOl-Shen 


io 




slla3o09_Shen 


A 




s20a3o05_LS 


cr 






s40a3o07_LS 


4 


4 


s40a3ol3_Shon 


4 


4 


Mur \7l\ 






20_3.8 


9 


4 


40_10.0 


4 


4 


40_13.0 


21 


15 


15_3.2 


10 


9 


15_4.0 


4 


4 


Ott |12| 






nomotol3 


143 


4 


nomotolS 


45 


4 


S15.0WHW02 


4 


4 


S20.0WHW02 


4 


4 


S25.0WHW02 


4 


4 



In Fig. [5j the waveforms associated with the acous- 
tic mechanism (Ott catalog) require the smallest SNR, 
followed by those of the neutrino mechanism (Mur cata- 
log) and the magnetorotational mechanism (Dim catalog). 
This hierarchy in minimum SNR, however, is not funda- 
mental, but a consequence of the fact that we have chosen 
to carry out this test using 7 PCs for each waveform cat- 
alog. Since the Ott catalog comprises only 7 waveforms, 
the set of 7 PCs completely spans it and allows perfect 
reconstruction, maximizing p(Z)|Ms) (Eq. [s]). In the case 
of less than perfect knowledge of the signal, the minimum 
SNR will always be greater. This is why the Dim and Mur 
catalogs, which have many more than 7 waveforms, re- 
quire larger minimum SNR than the Ott waveforms. 

The situation is somewhat different, when we re- 
calculate log Bsn using only 3 PCs. The Ott catalog is 
very small and rather diverse in the time domain. Its first 
few PCs do not efficiently span the catalog and, when 
only the first 3 PCs are used, the minimum SNRs for 
waveforms poorly reconstructed with these PCs increases 
dramatically, as shown for the Ott waveforms nomotol3 
and nomotolS in Tab.lH Some of the Mur waveforms also 
exhibit increased minimum SNRs, indicating that there is 
significant time-domain variation that is not captured by 
the first 3 PCs. The large Dim catalog, on the other hand, 
is very efficiently spanned already by its first few PCs, 
as previously pointed out by |36l [39| , and the minimum 
SNR for waveforms from this catalog remains practically 
unchanged when going from 7 to 3 PCs in our analysis. 



C. Determining the Core-Collapse Supernova 
Explosion Mechanism 

The basic assumption of this study is that the neutrino, 
magnetorotational, and acoustic core-collapse supernova 
explosion mechanism have robustly distinct GW signa- 
tures. In this section, we test this assumption by inject- 
ing waveforms into simulated noise and running SMEE 
on the data using PCs of waveform catalogs represen- 
tative of the neutrino, magnetorotational, and acoustic 
mechanisms. If our assumption is correct and the GW 
signatures of these mechanisms are truly distinct, then 
SMEE should (i) yield the largest value of log Bsn when 
the set of PCs is used that corresponds to the mechanism 
the waveform is representative of, and, (ii) lo gBjj ( Eq.[3t 



should be positive (and larger than ~5; see j IV A| if tT 
injected waveform is most consistent with mechanism i, 
negative if it is most consistent with mechanism j, and 
near zero if the result is inconclusive. 

We carry out our SMEE calculations for events located 
at 0.2 kpc, 2 kpc, and 10 kpc and with 3 and 7 PCs. Betel- 
geuse {a Orionis), a red supergiant star of 15 — 20 Mq 
pl6] . is located 197 ± 45 pc from Earth [113 and will 
eventually explode as a Type II-P core-collapse super- 
nova. Hence, studying a potential event at 0.2 kpc will 
tell us what we can learn from the observation of GWs 
from Betelgeuse's collapse and supernova. 2 kpc is still 
nearby on the galactic scale, but the Galactic volume out 
to this radius already contains hundreds of supergiants, 
one of which may make the next galactic supernova |118j . 
Finally, 10 kpc is the fiducial Galactic distance scale and 
we consider it to state what could be inferred throughout 
the Milky Way. As in Sec. |IVB we carry out SMEE runs 
with 3 and 7 PCs to study the sensitivity of the results 
on the amount of knowledge about the injected waveform 
we grant SMEE. 

Table |ll] summarizes the results from log Bsn calcula- 
tions for five representative waveforms from the Mur (neu- 
trino mechanism). Dim (magnetorotational mechanism), 
and Ott (acoustic mechanism) catalogs. The larger the 
value of log Bsn, the greater the confidence that the data 
contain a signal consistent with the employed set of PCs. 
log Bsn < -47 (for the 7-PC case; -21 in the 3-PC case) 
indicates that the signal is more consistent with noise. 
The results show that it is indeed possible to clearly as- 
sociate any injected waveform with its catalog and, thus, 
select the explosion mechanism for a core-collapse super- 
nova out to at least 2 kpc. However, at the galactic scale 
(10 kpc), a significant fraction of waveforms representa- 
tive of the neutrino mechanism, due to their intrinsically 
low GW amplitudes, can be told apart from noise only 
marginally or are even most consistent with noise. Typ- 
ical magnetorotational explosions and explosions driven 
by the acoustic mechanism are still clearly identifiable at 
10 kpc. However, when the reduced set of 3 PCs is used, 
the Ott catalog, which is not well spanned by only 3 PCs, 
suffers most and two out of the five representative Ott 
waveforms are not or only marginally identifiable with 
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TABLE II: logBsiv with respect to the Dim, Mur, and Ott PCs computed for representative injected waveforms from the 
Dim, Mur, and Ott catalogs, which we take to be representative of the magnetorotational, neutrino, and acoustic explosion 
mechanism, respectively. Results for source distances of 0.2, 2, and lOkpc are given. Each table entry shows the results for 3 
PCs (to the left of the vertical divider |) and the 7-PC result (to the right of the divider), log Bsn < —47 when 7 PCs are used 
and logi?5jv < —21 when 3 PCs are used indicates that the injected signal is likely consistent with noise while larger values 
suggests that the signal belongs to the signal model whose PCs were used in the analysis. 



Waveform 




log Bsn 




log Bsn 










log Bsn 






Name 




Dim PCs 




Mur PCs 










Ott PCs 








0.2 kpc 


2 kpc 


10 kpc 


0.2 kpc 2 kpc 


10 kpc 


0.2 kpc 


2 kpc 


10 kpc 


Dim ilTJ 
sllalo01_Shen 


38982|48686 


364|430 


-101-34 


9178125647 651203 


-22 


-42 


-25 


222 


-271-51 


-26 


-52 


slla3o09_Shen 


5x10'^ SxlO'' 


52403 54505 


2071 2129 


36170 58531 335 532 


-11 


-29 


23 


9191 


-26 39 


-26 


-49 


s20a3o05_LS 


3x10** 3x10" 


25738|26859 1005|1023 


4139 29285 15 240 


-24 


-40 


1 


2242 


-26-31 


-26 


-51 


s40a3o07_LS 


IxlO'^ IxlO'^ 


1x10^ 1x10^ 4737 5411 


47029 1x10^ 444951 


-7 


-12 


519 


56978 


-21 517 


-26 


-29 


s40a3ol3_Shen 


2x10^^ 2x10* 


2x10^ 2x10^ 8174 8302 


2x10-'^ 4x10^ 1830 3968 


48 


109 


160 


19781 


-25 145 


-26 


-44 


-— : Ts^n 1 

Mur |71| 
























20_3.8 


42051 10631 


16150 


-241-50 


5926414x10^ 56613775 


-2 


101 


-22 


-5 


-271-53 


-26 


-52 


40-10.0 


280 2362 


-23-33 


-26-53 


4x10^15x10^^ 419814981143 


149 


-26 


-29 


-27-54 


-26 


-52 


40_13.0 
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-26-54 
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-6 20 
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-26 
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Ott 12 
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-21111893 
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-25-53 
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295 65635 


-13 
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578 2417 
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-26-53 
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-25 
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-33 


-25-51 


1918 2969 -7-24 


-25 


-51 


1x10* 


1x10* 


1x10® 1x10® 


45853 


45829 



the 3 first Ott PCs at 10 kpc. 

Provided that there is confidence that a signal has been 
detected, we can compute log Bij (Eq. Isl) to study if the 
signal is more likely to be consistent with mechanism i or 
mechanism j. Since we know log Bsn for all signal mod- 
els, we can simply compute log Bij — log Bi^ — log Bj^ 
from the data for representative waveforms provided in 
Tab. HH We use 7 PCs for these calculations. 

In Fig. |6j we show results of injection studies of all 
waveforms from the Dim, Mur, and Ott catalogs run 
through SMEE and analyzed with the Dim, Mur, and Ott 
PCs at a source distance of 10 kpc. The full numerical re- 
sults on whose basis Fig.[6]was generated are available on- 
line |119j . The left panel depicts the log BDimMur result for 
injected waveforms from the Dim and Mur catalogs, that 
we take to be representative of the magnetorotational 
and neutrino mechanism, respectively. Even at 10 kpc 
the vast majority of waveforms characteristic for magne- 
torotational explosions are clearly identified as belonging 
to this mechanism. For the neutrino mechanism, the evi- 
dence is generally significantly weaker and only ~44% of 
the Mur waveforms are identified with logBDimMur < —100 
and none have logBoimMur < —1000, while ~19% are in 
the inconclusive regime of — 5 < log i?DiiiiMur < 5. 

In the center panel of Fig. |6j we show results for 
log i?Dimott for injected waveforms corresponding to the 
magnetorotational (Dim) and the acoustic (Ott) mech- 
anism. The case is clear cut and most waveforms are 



correctly identified as most likely belonging to their re- 
spective catalog/mechanism. Finally, the right panel of 
Fig. [6] presents log i?Murott for waveforms representative 
of the neutrino (Mur) and acoustic (Ott) mechanism. As 
in the previous panel, SMEE associates the waveforms 
corresponding to the acoustic mechanism with high con- 
fidence to the Ott catalog. The evidence suggesting cor- 
rect association of the neutrino mechanism waveforms is 
considerably less strong, but log i^Murott is still conclusive 
for ~88% of the Mur waveforms. 

Figure [7] shows the results for logSDimMur, fogSDimott, 
and log -BMurott obtained by SMEE with 7 PCs at a source 
distance of 2 kpc. Here, all acoustic mechanism wave- 
forms (Ott catalog), all magnetorotational mechanism 
waveforms (Dim call), and all neutrino mechanism wave- 
forms (Mur catalog) are correctly identified as belonging 
to their respective catalog and explosion mechanism. 



D. Deciding Between Rotating Accretion-Induced 
Collapse and Rotating Iron Core Collapse 

The waveforms of the Dim catalog are representa- 
tive of the GW signal emitted by rotating collapse and 
bounce of iron cores of massive stars with ZAMS masses 
> 8 — 10 M0. In the accretion-induced collapse (AIC) 
of rapidly rotating 0-Ne white dwarfs, very similar dy- 
namics occurs and the corresponding GW signals, as pre- 
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FIG. 6: Histograms describing the outcome of signal model comparisons by means of the Bayes factors logB^j = logp(D|A/i) — 
\o%p(D\Mj), where i i= j and Mi and Mj are signal models described by the Dim (magnetorotational mechanism), Mur (neutrino 
mechanism), and Ott (acoustic mechanism) waveform catalogs. The Bayes factors are computed with 7 PCs and for a source 
distance of lOkpc. A positive value logB^j indicates that the injected waveform most likely belongs to model Mi, while a 
negative value suggest that model Mj is the more probable explanation. The bars are color-coded according to the type of 
injected waveform. The results are binned into ranges of varying size from < —10000 to > 10000 and the height of the bars 
indicates what fraction of the waveforms of a given catalog falls into a given bin of log Bij. We consider the range of (—5, 5) of 
logSij as inconclusive evidence (see SIVA I. 
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FIG. 7: Same as Fig. [6] but computed for a source distance of 2kpc. 



dieted by Abdikamalov et al. [46] , share many of the ba- 
sic features of the rotating iron core collapse and bounce 
waveforms of, e.g., the Dim catalog (see the discussion in 
Sec. IV. C. of [13). Hence, it is interesting to see if our 
SMEE model selection algorithm can tell them apart. 

We compute the PCs for the Abd catalog in the same 
fashion as done previously for the Dim, Mur, and Ott cata- 
logs and inject all Abd and Dim waveforms into simulated 
Advanced LIGO noise. SMEE is then run with 7 PCs 
to calculate logSAbdDim- The result is shown in Fig. [8] 
for source distances of lOkpc and 2kpc. Full numerical 
results are available on-line |119| . 

In spite of the strong general similarity of rotating iron 
core collapse and rotating AIC waveforms, SMEE cor- 
rectly identifies the vast majority of injected waveforms 
as most likely being emitted by a rotating iron core col- 
lapse or by rotating AIC. However, for a source at lOkpc 
(left panel of Fig.[^, -6% of the Dim and -5% of the Abd 



are incorrectly identified as belonging to the respective 
other catalog. For an additional 2% of the Dim waveforms 
and 14% of the Abd waveforms, the evidence is inconclu- 
sive. 

At a source distance of 2 kpc (right panel of Fig. |8| , 
88% of the AIC (Abd) and 93% of the rotating core col- 
lapse (Dim) waveforms are correctly identified. 

If one placed trust in the reliability of less dominant 
and more particular features of waveforms in the underly- 
ing catalogs, one could use a larger number of PCs in the 
analysis. In order to study the effect of using an increased 
number of PCs, we re-run the Abd vs. Dim comparison 
with 14 PCs and find that the result is significantly worse 
than with 7 PCs: ~61% of the Abd waveforms and ^^23% 
of the Dim catalog are now incorrectly attributed to the 
respective other catalog at 10 kpc. This counter intuitive 
and at first surprising result is readily explained by the 
overall great similarity of the AIC and iron core collapse 
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FIG. 8; Outcome of the SMEE analysis of injected rotating iron core collapse (Dim catalog) and rotating accretion-induced 
collapse (AIC, Abd catalog) waveforms. The left panel shows results for a source distance of lOkpc and the right panel depicts 
the results for a distance of 2kpc. The Bayes factors logSAbdoim are computed on the basis of 7 PCs from the Abd and Dim 
catalog. A positive value of log BAbdDim indicates that an injected waveform is most likely associated with rotating AIC and a 
negative value suggests it to be more consistent with rotating iron core collapse. The results are binned into ranges of varying 
size from < —10000 to > 10000 and the height of the color-coded bars indicates what fraction of the waveforms of a given 
catalog falls into a given bin of log BAbdoim. We consider the range of (—5, 5) of log_Bij as inconclusive evidence (see SIVA I. 



waveforms and the nature of PCA and SMEE's Bayesian 
model selection. The most robust features of each wave- 
form catalog are encapsulated in its first few PCs. The 
first Dim and Abd PCs are indeed significantly different, 
but subsequent Abd and Dim PCs exhibit rather simi- 
lar secondary features. Since each PC carries the same 
weight in SMEE's evidence calculation, including a larger 
number of PCs dilutes SMEE's judgment in this case and 
leads to the observed false identifications. 



E. Testing Robustness with Unknown Waveforms 

In the previous sections, we have demonstrated 
SMEE's ability to identify an injected trial GW signal as 
belonging to a particular physical model (i.e., emission 
mechanism and/or explosion mechanism). In this, how- 
ever, we have drawn the injected waveforms directly from 
the catalogs used to generate the PCs. In other words, 
we have given SMEE (limited) advance knowledge about 
the injected waveforms. 

Here, we carry out a much more stringent test of 
SMEE's ability to select between models of the core- 
collapse supernova mechanism by injecting waveforms 
that were not employed in the initial PC generation 
and/or stem from completely independent catalogs. 



(see Sec. II Bl). The results of the log Bsn calculation 



for the magnetorotational, neutrino, and acoustic mecha- 
nism signal models are summarized in Tab. Ill DimExtra 



waveforms are identified as being most consistent with 
the Dim catalog and, hence, the magnetorotational mech- 
anism. This is true with high confidence when 3 or 7 PCs 
are used and for all DimExtra signals out to distances 
>10kpc. 

The Sch waveforms were generated with a completely 
different numerical code and thus allow for a truly inde- 
pendent test of our approach in SMEE. Also, unlike the 
Dim waveforms, the Sch waveforms are based on 3D sim- 
ulations. Hence, they are not linearly polarized. For con- 
sistency with our current approach, we neglect hx and 
inject only /i+ as seen by an equatorial observer. Re- 
sults of SMEE log BsN calculations for all injected Sch 
waveforms are summarized in Tab. |III[ SMEE correctly 
identifies all injected Sch waveforms as indicative of mag- 
netorotational explosions at a source distance of 2 kpc. 
At 10 kpc, still 91% of the injected Sch waveforms are 
attributed to the magnetorotational mechanism, which 
is an indication of the robustness of the GW signal asso- 
ciated with rapid rotation and magnetorotational explo- 
sions. The very few Sch waveforms that SMEE is not able 
to clearly associated with the magnetorotational mecha- 
nism have such weak SNRs that they are more consistent 
with noise than with any of the catalogs at 10 kpc. 



1. Magnetorotational Mechanism 

For the magnetorotational mechanism, we employ 
three additional Dim waveforms (DimExtra, Sec. IIBl) 
that were not included in the calculation of the Dim PCs. 
Furthermore, we inject waveforms from rotating mod- 
els of the Sch catalog of Scheidegger et al. jl(M 1107] 



2. Acoustic Mechanism 

We test SMEE's ability to identify core-collapse super- 
novae exploding via the acoustic mechanism by injecting 
the three DttExtra waveforms (see Sec. IICll. The re- 
sults of this test are again summarized in Tab. |III[ They 
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TABLE III: log Bsn for gravitational waveforms that were not included in the catalogs used for PC computation. The DimExtra, 
Sch, OttExtra, and Yak waveforms are discussed in ijll] The MurRem waveforms are three randomly selected waveforms from the 
Mur catalog that were removed before re-computation of the Mur PCs. Results are shown for source distances of 0.2 kpc, 2kpc, 
and 10 kpc and for evaluations using 3 PCs (to the left of the vertical divider |) and 7 PCs (to the right of the divider). Larger 
values indicate stronger evidence that the waveform is matched to the model catalog from which the PCs were constructed. 
logSsjv < —47 when 7 PCs are used and logBgiv < —21 when 3 PCs are used indicates that the injected signal is likely 
consistent with noise while larger values suggests that the signal belongs to the signal model whose PCs were used in the 
analysis. 
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suggest that the a-priori unknown OttExtra waveforms 
can be identified as belonging to the acoustic mechanism 
out to 2 kpc with great confidence when 7 PCs are used in 
the analysis. At 10 kpc, the waveforms are still correctly 
attributed to the acoustic mechanism, but the evidence 
is much weaker in the 7-PC case while the waveforms 
are more consistent with noise when the analysis is per- 
formed with only 3 PCs. The OttExtra 3 waveform, 
which is clearly identified at 10 kpc, has an extreme SNR 



of ^2530 at this distance, while the two other waveforms 
have SNRs of --50. SMEE's difficulty is illustrated in the 
right panel of Fig. |9j which indicates that the OttExtra 
waveforms reach the threshold of log Bsn > —47 only 
for SNRs >35, whereas Ott waveforms are identified al- 
ready at SNRs >4, if the full set of 7 PCs is used. This 
is a strong indication that the range of possible wave- 
form features associated with the acoustic mechanism is 
not efficiently covered by the 7 PCs generated from the 
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FIG. 9; Mean and 1-a range of logSsjv as a function of signal-to-noise ratio SNR comparing signal with noise evidence. 
The horizontal lines mark the threshold values of log Bsn above which we consider an injected waveform to be distinct from 
Gaussian noise. Left panel: Results for the Sch and DimExtra. These two were both reconstructed with 7 Dim PCs. Right 
panel: Results for the Yak, MurRem and OttExtra waveforms as reconstructed with 7 Mur for the first two and 7 Ott PCs for 
the latter. The Dim PCs very efficiently reconstruct the Sch and DimExtra waveforms at moderate SNRs while the Yak and 
OttExtra require very high SNRs to be distinguished from noise by the Mur and Ott PCs, respectively. 



Ott catalog. This could simply be attributed to the very 
small number of waveforms in this catalog. However, 
when studying the Ott and OttExtra waveforms, one 
immediately notes that the time between the first peak 
(associated with core bounce) and the second peak (the 
global maximum, associated with the non-linear phase 
of the protoneutron star pulsations) varies significantly 
between waveforms. Since we compute PCs in the time 
domain, such large-scale features are imprinted onto the 
PCs and make it difficult to identify waveforms whose 
two peaks are separated by significantly different inter- 
vals. An alternative method that may work much better 
for waveforms of this kind is to compute PCs based on 
waveform power spectra, which would remove any poten- 
tially problematic phase information. 



3. Neutrino Mechanism 

We test SMEE's ability to identify GW signals emit- 
ted by core-collapse supernovae exploding via the neu- 
trino mechanism in two ways. First, we remove three 
randomly selected Mur waveforms (Mur waveforms 20_3.4, 
12_3.4, and 15_3.2, labeling this set as MurRem) from the 
Mur catalog, recompute the PCs without these 3 wave- 
forms, then run SMEE to compute log Bsn for the three 
MurRem waveforms. The results, listed in Tab. IIIIl show 



at 10 kpc two out of three MurRem waveforms are picked 
out of the noise (though with relatively weak evidence) . 
This is consistent with the overall results for waveforms 
belonging to the Mur catalog discussed in Sec. |IVC| 

However, a comparison of the right panel of Fig. |9]with 
Fig. [5] shows that the MurRem waveforms require an SNR 
that is more than twice as high to reach values of log Bsn 
at which we can consider them to be distinct from Gaus- 
sian noise. This is most likely due to the rather large 
diversity of Mur waveforms. Components of relevance to 
the MurRem waveforms are apparently not captured in the 
first 7 PCs when these waveforms are not included in the 
PCA. 

A yet more stringent test is enabled by the waveforms 
of the Yelk catalog (see Sec. II A 1 1 that were obtained 
with a completely different numerical code. We inject 
the three available Yak waveforms into Advanced LIGO 
noise and run SMEE on them to compute log Bsn- We 



that SMEE is able to correctly identify MurRem wave- 
forms as GW signals consistent with the Mur catalog with 
strong evidence out to a distance of ~2 kpc and even 



list the results in Tab. HI SMEE correctly and clearly as- 
sociates the Yak waveforms with the Mur PCs at 0.2 kpc. 
At 2 kpc, the association is still possible, but at 10 kpc the 
Yak waveforms appear to be most consistent with noise 
for SMEE. The right panel of Fig. [9] shows that the Yak 
waveforms require an SNR to be clearly associated with 
the neutrino mechanism that is >7 times higher than for 
MurRem waveforms and more than ~17 times higher than 
for Mur waveforms. This rather disappointing result can 
be explained as follows: While the Yak waveforms are 
qualitatively very similar to the Mur waveforms, they dif- 
fer significantly in quantitative aspects. The Yak wave- 
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forms are generally only half as long ('^Is for Mur and 
0.5 s for Yak, whose models explode much earlier than 
the Mur models). Furthermore, the Yak waveforms have 
considerably more power at frequencies above ^800 Hz 
and their energy spectra peak at ^1000 Hz while most 
of the emission in the Mur waveforms occurs at or below 
~400Hz. This may be due to the more simplified treat- 
ment of gravity and neutrino microphysics and transport 
in the study of Murphy et al. [7T] underlying the Mur 
catalog compared to the work of Yakunin et al. [58] that 
led to the Yak catalog. 



V. SUMMARY AND DISCUSSION 

In this article, we have described the Supernova Model 
Evidence Extractor (SMEE), a novel Bayesian approach 
to inferring physical information from observations of 
GW bursts emitted in stellar collapse and core-collapse 
supernovae. SMEE decomposes catalogs of simulated 
GW signals into their principle components and employs 
the Nested Sampling algorithm to compute the evidence 
that a trial signal injected into GW detector noise be- 
longs to a given catalog. 

It is evident that core-collapse supernovae are pow- 
ered by the release of gravitational energy in gravita- 
tional collapse, but the central unsolved problem of core- 
collapse supernova theory is by what mechanism this en- 
ergy is transferred from the collapsed core to the stel- 
lar mantle to drive the explosion. Simulations show that 
multi-dimensional GW-emitting dynamics is a crucial in- 
gredient to all potential explosion mechanisms. Hence, 
GWs could be ideal probes for the core-collapse su- 
pernova mechanism, provided that different mechanism 
models have clearly distinct GW signatures. In this pa- 
per, we have considered (i) the neutrino mechanism, {ii) 
the magnetorotational mechanism, and {in) the acous- 
tic mechanism for core-collapse supernova explosions and 
have treated mechanism as equally probable (i.e., having 
the same prior probability). The primary and dominant 
multi-dimensional dynamics and GW emission processes 
of these mechanisms are, ordered in the above order 
of mechanisms, (i) convection/turbulence and accretion 
downstreams, {ii) rapid rotation, and {Hi) protoneutron 
star core pulsations, respectively. 

Using GW signal catalogs based on simulations rep- 
resentative of these three mechanisms, we have applied 
SMEE to infer the explosion mechanism underlying trial 
waveforms injected into simulated Advanced LIGO noise. 
Our results show that our Bayesian approach is capable of 
identifying any of the considered explosion mechanisms 
with high confidence for core collapse events occurring 
out to a distance of ~2 kpc even with only rudimentary 
knowledge of the precise shape of the expected signal. 

Magnetorotational explosions, leading to particularly 
strong GW signals with very robust common features, 
can be clearly identified throughout the Milky Way 
{D > 10 kpc). Our results also suggest that it is 



possible to further distinguish between rapidly rotating 
accretion-induced collapse of massive white dwarfs and 
rapidly rotating iron core collapse for events occurring 
in the galaxy, provided that the differences between the 
GW signals predicted by current simulations of rotating 
accretion-induced collapse and rotating iron core collapse 
are reliable. 

GW signals emitted by neutrino-driven explosions have 
systematically lower amplitudes and, hence, are harder to 
distinguish from detector noise. Moreover, GW emission 
from convection/turbulence and accretion downstreams 
has not been as extensively studied, the available number 
of model waveforms is an order of magnitude smaller, and 
the currently predicted GW signals are not as reliable as 
in the case of rapidly rotating collapse. This reduces the 
efficacy of the principal component decomposition and, 
combined with the overall weakness of the signals, lim- 
its SMEE's robust reach for the neutrino mechanism to 
<2 kpc for the currently available set of model waveforms 
and Advanced LIGO in broadband configuration. 

The GW signals from core-collapse supernovae driven 
by the acoustic mechanism are strong and would likely 
be detectable by Advanced LIGO throughout the Milky 
Way. However, the set of available model waveforms for 
this signal type is very limited and individual waveforms 
differ significantly at large scales in the time domain 
while having similar frequency content. The inclusion 
of phase information (by computing PCs in the time do- 
main and using their full complex Fourier transforms) in 
the current incarnation of SMEE is sub-optimal for such 
signals and the PCs only inefficiently span the space of 
possible waveforms. This makes it difficult for SMEE to 
clearly identify the acoustic mechanism for core collapse 
events occurring at distances significantly greater than 
-2 kpc. 

This study is the first systematic attempt at inferring 
core-collapse supernova physics from observations of GW 
bursts emitted by multiple different underlying mecha- 
nisms. While this is a significant step beyond previous 
work that focused only on GWs from rotating core col- 
lapse [32l [39] , our present study still suffers from a num- 
ber of simplifying assumptions: We have considered only 
a single detector, Gaussian noise, and linearly polarized, 
optimally oriented GW emission. Real core-collapse su- 
pernovae will emit in both GW polarizations and will 
have arbitrary, generally non-optimal orientation with re- 
spect to observatories on Earth. Advanced LIGO class 
GW detectors will operate as networks and observations 
will be coincident between three or even four detectors, 
and the noise backgrounds of these detectors will not be 
stationary and Gaussian, but non-stationary and glitchy. 
We have also been agnostic with regard to the prior prob- 
ability of the three considered core-collapse supernova 
mechanisms. Input from astronomical observations and 
theory could be used to generate (approximate) prior 
probabilities for each mechanism, which could then be 
included in a future analysis. 

Near future work on SMEE will be directed towards 
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incorporating detector networks, variations in detector 
configurations and sensitivities, both GW polarizations, 
non-Gaussian noise, arbitrary source-detector orienta- 
tions, and improved principal component analysis with 
and without reliance on the signal phase. However, even 
with these improvements, successful extraction of physics 
from core-collapse supernova GW signals will crucially 
depend on the availability of extensive catalogs of reliable 
predictions for both h+ and hx ■ These must be provided 
by the core-collapse supernova modeling community on 
the basis of full 3D simulations. 
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